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We address the low-energy effective Hamiltonian of electron doped d° perovskite semiconductors 
in cubic and tetragonal phases using the k • p method. The Hamiltonian depends on the spin-orbit 
interaction strength, on the temperature-dependent tetragonal distortion, and on a set of effective- 
mass parameters whose number is determined by the symmetry of the crystal. We explain how 
these parameters can be extracted from angle resolved photo-emission, Raman spectroscopy, and 
magneto-transport measurements and estimate their values in SrTiOs. 
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I. INTRODUCTION 

Transition metal oxides with perovskite structures 
exhibit a wide variety of interesting and often useful 
effects including colossal magnetoresistanceP high T c 
superconductivity,^ and ferroelectricit}^. Correspond- 
ingly, these materials have received intense experimental 
and theoretical attention for over half a centur^. Within 
the perovskite family, the d° materials have received par- 
ticular attention, often because of their large band gaps. 
SrTiOs, for example, is perhaps the most common sub- 
strate for the epitaxial growth of oxide materials. Re- 
cently there has been growing interest in the transport 
properties of lightly electron doped d° perovskites P In 
KTaOs, for example, strong spin-orbit (SO) coupling 
facilitates electrical manipulation of spin in a field ef- 
fect transistor geometry. 6 The two-dimensional electron 
systems which form at interfaces between d° materials 7 
show intriguing magnetic phases 8 and peculiar magneto- 
transport features Advanced epitaxial growth tech- 
niques enable J-dopingof oxides^ and the fabrication of 
oxide hetrostructuresP^ These relatively recent rapid ad- 
vances could, it is hoped, eventually lead to useful oxide 
based nano-electronic devices 

The low-energy band structure of an oxide provides a 
starting point for understanding not only its bulk trans- 
port characteristics but also its electronic properties near 
£-doped layers and near interfaces. First principles elec- 
tronic structure theory methods^^ are usually efficient 
for determining the gross structure of a band, but are 
not sufficiently accurate to nail down the fine features 
that determine the electronic properties of the states 
at the bottom of the conduction band that are impor- 
tant in weakly doped bulk materials, and in low-carrier- 
density two-dimensional electron systems. In particu- 
lar, it appears that at present bulk band structures in 
d° perovskites are not known accurately enough to pre- 
dict the two-dimensional bands of J-doped oxides or of 
interface-localized bands in oxide based het restructures. 
This paper is primarily motivated by the goal of assisting 
progress in this direction. 

The k • p methocP^El] offers an alternative and a po- 
tentially more accurate route for characterizing band 



structure near the conduction band minimum. The 
method provides an effective Hamiltonian that depends 
on a set of phenomenological parameters which can be 
small in number when band extrema occur at high- 
symmetry points in momentum space. The utility of 
this method hinges on the ability to extract accurate 
parameter values from experiments. In the case of per- 
ovskites the most valuable experimental probes are angle 
resolved photo-emission (ARPES), Raman spectroscopy, 
and magneto-transport measurements. 

Many of the most studied oxides have conduction-band 
minima located at the center of the Brillioun zone. We 
therefore apply the k • p method to obtain an effective 
low energy Hamiltonian near the T point. At high tem- 
peratures, perovskites typically have cubic symmetry. 
As the temperature is decreased the symmetry is usu- 
ally lowered, most commonly to either orthorhombic or 
tetragonal. The distortion can be driven by the motion 
of atoms along one of the cubic axes (e.g. in BaTiOs) or 
by a rotation of the oxygen octahedras (e.g. in SrTiOs). 
Structural phase transitions can also be induced by ap- 
plied stress p3 

In this work we focus on the cubic and tetragonal 
phases. In section |H] we briefly describe the k • p method 
and then use it to derive the low energy effective theory 
of a dP perovskite in the vicinity of the Y point. In section 
III we elaborate on experimental methods for obtaining 
the parameters of the k • p Hamiltonian. Using the ex- 
perimental data accumulated over the past few decades 
we then study the effective Hamiltonian of the conduc- 
tion bands of SrTiOs in Section |IV[ We summarize in 
section [V| 



II. LOW ENERGY THEORY 

For many perovskites of current interest such as 
SrTiOs the conduction band minima is at the Brillouin- 
zone center T-point. For momenta near the T-point the 
crystal field splits the ten <i-bands into four high energy 
e g bands, and six lower energy t2 g bands. Because the 
crystal field induced gap is typically a few eV's, it is suf- 
ficient to consider the t2 g bands when constructing a low 
energy theory of weakly-doped d° materials. In the cu- 
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bic phase the t<i g bands are degenerate at the T-point if 
spin-orbit interactions are neglected, but are weakly-split 
by typical tetragonal or orthorhombic distortions and by 
weak spin-orbit interactions. Unless the Fermi energy is 
large compared to these splittings, spin-orbit and distor- 
tion related band parameters must be accurately known 
in order to achieve a reasonable description of electronic 
properties. 



perovskite structure by inversion symmetry. The matri- 
ces ft L and h so account phenomenologically for tetrago- 
nal distortion and SO interactions at the T point and are 
discussed more explicitly below. 

The wave functions at the zone center have no covalent 
character and can be spanned by the t^g basis 



(5) 



A. Effective Hamiltonian 

The unperturbed Hamiltonian in the k • p perturba- 
tion theor}ff^ is 
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Hq consists of three terms: the kinetic energy term, the 
lattice potential term V(r), and the spin-orbit term (a is 
the Pauli matrix vector). The k • p Hamiltonian, which 
acts on the periodic part of the Bloch state, includes a 
second term which accounts for the dependence of band 
wavefunctions on Bloch wave vector k: 
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The k • p method exploits the high symmetry at the T 
point to classify the k = wave functions by irreducible 
representations (irreps) of the appropriate point group 
symmetry. It then uses perturbation theory 
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to evaluate t<i g projected Hamiltonian corrections to sec- 
ond order in the Bloch wavevector k. Hereafter we use 
units in which h = 2m = 1 where rn is the bare mass of 
the electron. The six ti g band energies e(k) then follow 
from the secular equation 



det[ft so + K + ft(k) - e(k)J] = 0. 



(4) 



In Eq.(|3| {l^j)} label a basis set for the t^ g bands and (j) a 
is summed over bands outside the ti g manifold. The first 
order term was omitted in Eq. ([3| since it vanishes for the 



Here X, Y and Z correspond respectively to the \yz), \xz) 
and \xy) t^g orbitals. Below we obtain the Hamiltonian 
matrix in this basis. 

The lattice term ft L is non zero in the tetragonal phase. 
If we choose a convenient zero of energy and set the z axis 
along the tetragonal axis then ft L has a single non-zero 
matrix element: 



(Za\V\Za) = A T , 



(6) 



where a accounts for the spin. The SO term in the Hamil- 
tonian is 



(h so ) 



so)ia,kf3 



(Z i a\k-v\Zkf3) = (Z i \AMk)-(a\a j \/3), (7) 



where A ex W x p and & is one of the orbital ba- 
sis functions. Because A transforms as a pseudovector, 
(£i|Aj|£fc) ex eijk where is the third rank antisymmet- 
ric tensor. For example, (X\A Z \X) and (X\A Z \Z) vanish 
under reflection off the x-z plane. Furthermore, since the 
matrix elements ^ must be imaginary 
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Strictly speaking, SO coupling is described by two pa- 
rameters in the tetragonal phase. However we neglect 
this small correction since it is of order of A T over the 
band gap compared to the spin-orbit coupling term we 
retain. 

The k-dependent part of the Hamiltonian h is obtained 
using Eq.([3|. We show in the appendix |A| that 



h 



h 



(9) 



with 
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-A/45 k, 
Af^kyk 
M 4 (k 2 x + k 2 y ) 
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In the tetragonal phase the h matrix depends on eight bic phase parameter values become independent of their 
real parameters (only A/45 may be complex). In the cu- subscript labels (e.g. £4 = £5 —> L) and h then de- 
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pends on only three parameters. The energy dispersion 
relations follow from Eqs.(4jG — 10). Because the Hamil- 
tonian is time-reversal invariant and has inversion sym- 
metry it gives rise to three doubly-degenerate bands. 

In the next section we discuss zone-center wave func- 
tions and energies. The wavefunctions play a crucial 
role in matrix-element considerations which powerfully 
expand the ability of ARPES experiments to determine 
the parameters of the k • p-Hamiltonian. The zone-center 
energies can be compared with t<i g band-splitting values 
obtained by Raman spectroscopy. 



B. Zone center energies and wave-functions 

The Hamiltonian at the zone center is h L + h so . The 
energies are therefore 
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(Energy has been shifted so that will vanish.) In the 
cubic phase the t<i g bands transform as in the absence 
of spin-orbit coupling. SO interactions split the bands to 
r++T+. When there is a tetragonal transition, the four- 
fold degenerate states further split to I\ + IV The 
notation in Eqs.(ll) correspond to these latter irreps. 



The (unnormalized) wave functions corresponding to 
the energies (11) are 
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where 



D = 3A x /2 - A so /2. 



(13) 



(14) 



It is interesting to follow the evolution of the bands as 
the ratio between A T and A so is varied from zero to 
infinity. The two limits are given in table [TJ In the cubic 
phase the states j^ 6 ,^ 70 } are degenerate and are spilt 
off from the {^ 7b } states by an energy of A so . In the 
tetragonal phase when |A so | > |A T | the states group to 
the three doubly degenerate pairs ^ 6 ,^ 7a and ip 7b . As 
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(x t + iY t ),r 6 (r 8 ) 


(Xi. - iYi.) ,r 6 (v 5 ) 
(X r + iY r ),r 6 (T 5 ) 


[X f -iY t + 2Z i ],T 7 {T s ) 

[x l + iY l -2Z t ],r 7 (r 8 ) 

[-X t + iY r + , r 7 
[Xi + iYi + Zf],^ 


z u r 7 (r 4 ) 
-z t ,r 7 (r 4 ) 
-(x t -iy t ),r 7 (r 5 ) 
+ in), r 7 (r 5 ) 



TABLE I: Zone center wave functions in the cubic phase with 
SO interactions (left column) and in the tetragonal phase in 
the absence of SO interactions (right column). 



the temperature is lowered the four ip 7 states mix. If 
eventually |A so | <C |A T | then the and ifj 7b states 
combine to give the Z± state which is purely tetragonal 
in character. 

In the following section we discuss energy dispersion re- 
lations along symmetry lines and planes, which can be di- 
rectly related to ARPES measurements and enable some 
qualitative insights into the relationships between Hamil- 
tonian parameters and the field-orientation dependence 
of magnetoresistance-oscillation frequencies. 



C. Energy dispersion relations for high- symmetry 
lines and planes 



In general Eq.Q must be diagonalized numerically. 
However, simple energy dispersion relations exist along 
high symmetry directions and in high- symmetry planes. 

When the tetragonal distortion is large and SO inter- 
actions can be neglected, the ti g bands split into T4 + r$ 
bands. In this limit (to order & 4 /A T ) 



e 4 (k) = A T + M4k\ 
e 5 ±(k) = B + k 

± 



2 ■ M±kl 



(15) 



-Mi] k%kl 



where k 2 = k 2 x + k 2 y and B± = (£ 5 ± M\)/2. To leading 
order in A S o, the energies remain unchanged whereas 
the r 5 ± energies vary linearly in opposite directions. The 
energies (15) are valid for any value of A T (but still ne- 
glecting A so ) in the k z = plane. Similarly for the 
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k y = plane 
e 5 _(k) = M\k 2 x +Mikl 

£4,5+ (k) = 



At Mi + £5 k 2 , U + Mi fc2 



± ^ [A T + (X 4 - £ 5 ) k 2 x + (£ 4 - A*£) k 2 z ] : 



+4|^45| 2 fc^ 2 
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The k • p Hamiltonian for the ti g bands in the cubic 
phase is identical to that of the valenc e band p-states of 
zinc-blende type semiconductors! 17 * 18 ^ In the presence of 
moderate SO interactions the dispersion relations along 
the three equivalent principle axes are 



e 7 (k) 
e s± (k) 



Mk 2 
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For strong SO interactions the ip 7 and ^ 8 states can be 
approximately decoupled to order k 4 /A so . The energy 
dispersions are then 

e 7 (k) = A so +Ak 2 (18) 
e 8 (k) = Ak 2 ± ^B 2 k± + C 2 (k 2 k 2 + k 2 x k 2 z + k 2 k 2 ) 

where A = 1 + (L + 2M)/3, B = (L - M)/3, and C 2 = 
[TV 2 — (L — M) 2 ] /3. Expressions (18) were obtained by 
Dresselhaus et. al.[T7]. 

ARPES measurements are frequently set to measure 
the energy dispersion in the k x — k y plane. For k z = 
the dependence of band energies on momenta is similar 
in the dominant tetragonal-splitting and dominant spin- 
orbit coupling limits (compare Eqs.(15) and (18)). One 
way to determine which of the two interactions is domi- 
nant is to probe the dispersion relation along z. A second 
way is to monitor the evolution of the bands as a func- 
tion of temperature. Additional methods are explained 
in section Hill below. 

The parameters of the effective Hamiltonian in the 
tetragonal phase are temperature dependent. As T is 
lowered the tetragonal distortion increases and the en- 
ergy bands change accordingly. For some crystals, such 
as SrTiOs, the deformation is well described by a sim- 
ple order parameter^. It is then possible to express the 
temperature dependence of the different Hamiltonian pa- 
rameters via a single temperature dependent order pa- 
rameter. 



III. EXPERIMENTAL METHODS FOR 
DETERMINING HAMILTONIAN PARAMETERS 

The utility of the k • p method depends on the ability 
to extract accurate values for the Hamiltonian parame- 



ters from experiments. ARPES, magneto-transport, and 
Raman spectroscopy measurements are three of the most 
useful experimental probes for band parameters. In this 
section we focus on the ways in which these techniques 
can be exploited for d° perovskites with an emphasis on 
experimental signatures of the tetragonal distortion. 



A. Raman spectroscopy 

Raman spectroscopy is routinely used to measure the 
spectra of solids^. For a low doped d° perovskite Raman 
spectra can determine the band gaps at the zone center. 



As explained in section |IIC| distinguishing between A T 
and A so using ARPES measurements may prove diffi- 
cult. The band gaps depend both on SO interactions 
and on the tetragonal distortion. Spectroscopically mon- 
itoring the energy gaps as a function of temperature and 
comparing with Eqs.(ll) provides in principle sufficient 
information to determine A so and A T . 



B. ARPES 

Angle Resolved Photoemission Spectroscopy (ARPES) 
has now been developed into a widely applicable ex- 
perimental tool for the measurement of bulk and sur- 
face electronic states P3 In a typical measurement inci- 
dent monochromatic radiation excites electrons in occu- 
pied crystal states and unbinds them from the crystal. 
In the sudden approximation electrons are promoted di- 
rectly from a crystal state to a vacuum plane wave state. 
In this approximation the intensity of the ARPES sig- 
nal associated with in-plane electron momentum kn and 
energy u is 



J(k||,cj) oc VIM 



p,nk 



|M„(k|| )W )/H. (19) 



Here the z-axis is set perpendicular to the sample's sur- 
face and we assume that the photon energy is calibrated 
to probe the k z = plane. A n is the electron spectral 
function of band n, / is the Fermi distribution function, 
and 



M p , nk oc (p|A.p|tf nk ) = A.p(p|tf nk ) 



(20) 



gives the probability amplitude for an electron in an ini- 
tial state \I/ n k|| to transition to a plane wave state p via 
a photon field A. The photo-emitted electrons are se- 
lectively collected according to their emission angle and 
energy. Therefore in a given measurement the outgo- 
ing momentum p in Eq.(19) is fixed by the position of 
the detector and by the energy of the incoming photon. 
The component of the momentum parallel to the surface 
must equal the momentum of the initial state to within 
a surface reciprocal lattice vector. 

In principle with sufficient ARPES data the occu- 
pied energy bands can be accurately mapped. The k • p 
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Hamiltonian parameters can then be determined using 
the dispersion relations in section II C In practice, how- 
ever, experimental limits on energy and momentum res- 
olution combined with the relatively large number of 
Hamiltonian parameters and the possibility of surface 
states that obscure bulk bands, often complicate com- 
parisons between theory and experiment. 

As we now explain, additional band structure infor- 
mation can sometimes be drawn from systematics in the 
dependence of the ARPES matrix elements on the surface 
reciprocal lattice vector added to the transverse momen- 
tum. Matrix elements contributions from particular ti g 
orbitals frequently vanish at particular reciprocal lattice 
vectors either because of symmetry considerations or be- 
cause of photon polarizations. By noticing the reciprocal 
lattice vectors at which the signal from a particular band 
is absent or very weak, it may be possible to identify 
the tig components which contribute dominantly to that 
band. This orbital information strongly constrains the 
band model. 

In the sudden approximation 

<P|*nk,|> = *k||+G||, P || 5Z4 n) ( k ll) 
3 

x J drdze- i[ ( G w r+PzZ )^(r,z). (21) 

Here Gm is the surface-plane projection of a reciprocal 
lattice vector, and £j are the ti g basis functions given 
by Eq.([5| for the conduction band initial wavefunction: 



^ ?1 



exp (iku 



^3 a 3 



f\j. The ^-function in Eq.Q 



reflects the conservation of the in-plane crystal momen- 
tum in the photon assisted scattering process of the elec- 
tron. 

We illustrate the usefulness of the matrix element effect 
by considering M for G = ([00] BZ) and for G along 
the x-axis ([10] BZ). In the first case M always vanishes 
since all £j's are odd with respect to reflection off either 
the z-x or the z-y plane. There is therefore no ARPES 
signal in the [00] BZ for ti g conduction band states. For 
the [10] BZ 



M [oi] 



oc 



drdze- i( - G ' x+p ' x ^Y(r,z). 



(22) 



Contributions from other ti g components of vanish 
because of their reflection symmetry in the x-z mirror 
plane (see Figfl]). Therefore only wave functions con- 
taining a Y orbit al will be detected in this case. Recent 
experiment J 22 * 23 ! on bulk SrTiOa find a single (doubly 
degenerate) band for k z = in the [01] BZ in the cu- 
bic as well as in the tetragonal phase (see section |lV| ) . 
The Hamiltonian described by Eqs.(|6]— [l0|) then indicates 
that the ATs and A so must be sufficiently small so that 
any hybridization between the ^-orbitals is negligible, 
(see Figj2|. 

As evident from the Hamiltonian (|3| and illustrated in 
FigJ3]the tig d-orbitals are hybridized by N. The influ- 
ence of N is most pronounced along the main diagonals. 
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FIG. 1: (Color online) Left: High symmetry points and lines 
in the Brillouin Zone (BZ) of a simple cubic lattice. Right: 
Experimental geometry used in simulations of ARPES data 
in the [10] BZ. The photon source (green), with polarization 
in the xz-plane (blue), excites an electron to a high energy 
state that is emitted towards the detector (red dashed). Al- 
though the experimental geometry is unchanged by a reflec- 
tion through the xz-plane, the yz and xy orbitals are odd 
under this operation. This leads to measurement of only the 
zx band in this measurement. 



A(Y) r A(Y) S r £ 



A(X) r A(X) A(X) F A(X) 

FIG. 2: (Color online) Simulated ARPES signal in the [10] BZ 
for a temperature of 120K, L/M = 1/8, A so = 0, and N = 0. 
The resolution for all ARPES simulations was set to lOmeV. 
Only a single elliptical FS cross-section is seen (bottom left). 
Energy distribution curves (EDCs) have been included along 
several high symmetry directions - in this case showing only 
a single band associated with the xz basis state. 



For example in the [110] direction N induces a momen- 
tum dependent gap of 2Nk 2 . 

Spin-orbit interactions will also mix the ti g d-orbitals 
however in contrast to the N ^ scenario they have no 
preferential direction. When the SO splitting is larger 
than the Fermi energy, the ARPES spectrum along the 
E direction is similar to the spectrum in the N ^ 0, 
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FIG. 3: (Color online) Comparison of d-orbital mixing. Sim- 
ulated ARPES signal for a temperature of 120K, L/M = 1/8, 
and A T = 0. (a) For N/M = 3 and A so = the hybridization 
is most pronounced along E but unseen in the EDCs along 
A(X), and A(y). (b) For N/M = and A so = 3e F one 
band has moved above the Fermi energy. The hybridization 
between the basis states is seen in the EDCs along all direc- 
tions. This experimental feature can be attributed to the lack 
of a preferential direction of the SO interaction. 



A so = case. However, unlike the N ^ case the 
photoemission spectrum is also altered along the x and 
y directions. This is evident along k x in the simulated 
ARPES data of Figj3J where the induced hybridization 
of the basis functions cause the previously dark band to 
become visible. 

The SO and tetragonal energies A so and A T determine 
the band splitting at the Y point (see Eq.([n])). In the cu- 
bic phase the value of A so can be extracted directly from 
an ARPES measurement in the [10] surface BZ (see Fig. 
^l). This simple picture is complicated in the tetragonal 
phase. The case where A so > e F and A T < e F is readily 
distinguished from the opposite limit by analysis of the 
dispersion along k x . As evident from Figs. [4}d and [4^ 
only in the case where A so > e^, is the dark weakly dis- 
persive band visible away from the T point. When both 
A so and A T are less than the Fermi energy, the Energy 
Distribution Curves (EDCs) at the F point can be used 
to distinguish between the two cases. This can be seen 
in Fig. [5) 



C. Magnetic Oscillations 

Magnetic oscillations in various physical properties 
such as the conductivity (Shubnikov - de Haas effect) 
and the magnetic susceptibility (de Haas - van Alfen ef- 
fect) provide invaluable information on the band struc- 
ture of solids^H^. The frequency of the oscillations F 
is related to the extremal cross-sectional area Ak of the 
Fermi surface in a plane perpendicular to the magnetic 
field through the Onsager relation F = 0o^./c/47r 2 . Here 
0o = hc/e is the magnetic flux quantum. Measuring F 
as a function of charge density, magnetic field orienta- 
tion, and temperature also makes it possible in principle 
to determine all the phenomenological Hamiltonian pa- 
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FIG. 4: (Color online) Simulated ARPES measurement for 
the [10] BZ for different values of A so and A T . L/M — 1/8, 
and N = in all figures, (a) T = WOK, A so = 0.5e F , A T = 0. 
As evident from all EDCs the SO splitting hybridizes the basis 
states. If the temperature is lowered, inducing a structural 
phase transition that is large, one band moves above the Fermi 
energy. This relatively weak hybridization is seen in the Fermi 
surface (FS) of (b) where T = 20K, A so = 0.5e F , A T = 3e F . 
In contrast, for T = 20K, A S o = 3e F , A T = 0.5e F the strong 
hybridization of the basis states leads to a more symmetric 
FS (c). This feature is also seen by comparing the EDCs of 
(b) and (c). 



rameters. 

In the naive picture of three ellipsoidal decoupled d- 
bands the cross sectional areas are simply given by el- 
lipses. However this oversimplified scenario breaks down 
for any realistic system due to the hybridization of the 
d-orbitals by iV, and by the SO interactions. Avoided 
crossings of the overlapping energy bands then result in 
more complicated energy surfaces. 

To illustrate the variety of possible shapes of electron 
pockets we consider a simple case with a small but finite 
band mixing (e.g. N > 0). A T /e F = 0.5. 

The cross sectional areas for three high symmetry di- 
rections of the magnetic field are depicted in the top 
row of Figj6] for the tetragonal phase. As A T in- 
creases the most energetic band is gradually depleted 
and the electronic charge is redistributed amongst the 
other two Fermi pockets. Eventually for A T /e F > 1 — 
min(L/M, M/L) there is no band crossing between the 
X7/-band and the other two bands. 

Avoided crossings in the cubic phase result in non- 
elliptical cross-sections as well. The extremal cross- 



7 



(a) 



(b) 



MY) T AfY) £ r 



Am r A(Y) £ r s 



V 


\\ II 

V 


o o 


v 


V 


A 

( 1 




> 








v 

f 


-e, g 





V 



itx) r itxi i(X) r a<x> 




FIG. 5: (Color online) When A S o and A T are both small, 
although all three bands can be seen directly in a measure- 
ment of the [10] BZ it can be unclear how to extract these 
parameters. This is seen by comparing (a) where A S o = .5e F . 
A T = 0.3e F with (b) where A S o = 0.3e F . A T = 0.5e F . A 
careful analysis of the EDC at the F point can distinguish 
between the two scenarios. This is shown in (c) where the F 
point EDC is shown for (a) in dashed-red and for (b) in blue. 
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FIG. 6: (Color online) Extremal cross sectional areas for mag- 
netic field oriented along [001] (left) [010] (center) and [100] 
(right) . The extremal orbits have been organized by size from 
largest to smallest and shown as dashed-blue, dot-das hed-red, 
and dotted-green, respectively. Top row (a-c) corresponds to 
A T = 0.5e F , middle row (d,e) corresponds to A S o = 0.5e F , 
and bottom row corresponds to N = 0.5M. 



sectional areas along high symmetry directions are de- 
picted in Figj6] for A so /e F = 0.5 (center row) and for 
N/M 0.5 (bottom row). 

Our discussion ignores the possibility of multiple do- 
mains in the distorted state, and neglects magnetic 
breakdown. The latter is likely present in magnetic os- 
cillation measurements on these materials because of the 
close approaches between extremal cross-sections^ be- 
longing to different bands. 



IV. SrTiOs 

Bulk STO is a band insulator with an energy gap of 3.2 
eV. By chemical substitution of the Ti or the Sr atoms or 
by introducing oxygen vacancies it is possible to electron 
dope the system with a high level of precision. STO has 
cubic symmetry at room temperature, however at 105K 
it undergoes a antiferrodistortive structural transition to 
a tetragonal phase. Below the critical temperature neigh- 
boring TiOq octahedras continuously rotate in opposite 
directions by an angle of up to a few degrees. 

Although STO has been studied for many years, there 
are only a few experimental results that can shed light 
on the structure of its conduction bands. We therefore 
resort to a 5-parameter model in which the Hamiltonian 
is parameterized by A so ,A T ,M, iV, and L, i.e. h is ap- 
proximated by its cubic phase form. 

The experiments that do exist appear to partially con- 
tradict one another. Based on Raman spectro scopy and 
Shubnikov de-Hass measurements Uwe et. a/! 24 * 30 * con- 
cluded that A so « 18meV and A T ~ 1.5meV. On the 
contrary, Chang et. al. 23 using ARPES do not observe 
a SO induced gap at the zone center and conclude that 
A T « -25meV. 

Supporting evidence for the smallness of A so is pro- 
vided by the matrix element effect. Experiments 2 ^ ob- 
serve only, what should be according to our matrix el- 
ement analysis, the X orbital in the [10] BZ and th e Y 
orbital in the [01] BZ. As explained in section IIIB 1 the 
lack of hybridization between t^ g orbitals implies that 
both A so and N are very small. Additional proof that 
N < M, L is provided the ARPES EDCs which reveal no 
special features of the energy along E. In addition, these 
curves yield values for the effective masses from which it 
follows that 



M « 0.84 , L ^ 0.14. 



(23) 



Raman spectroscopy measurements^ find energy gaps 
of approximately 2meV and 18 meV between conduc- 
tion bands at the T point suggesting that A T and A so 
have very different magnitudes. The larger of the two 
scales can be identified as tetragonal or spin-orbit from 
the dependence of magnetic oscillation frequency F and 
cyclotron mass m* on density and field orientation. Fig(7| 
depicts the dependence of F and ra* on density and on 
A T . The dependence can be expressed through a single 
parameter A T /e F if F is scaled with n 2//3 where n is the 
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FIG. 7: (Color online) Dependence of magnetic oscillation 
frequency and cyclotron mass on A T . Here we set L/M — 1/8, 
N = 0, and A S o = 0. (a,b,c): Scaled SdH frequency for 
magnetic fields along [001], [110], and [100] as a function of 
At. (d,e,f): Cyclotron mass as a function of A T . 



In this work we used the k • p theory to construct the 
general low energy theory for the conduction bands of 
these materials both in the cubic and in the tetragonal 
phases. We then employed the theory to estimate the 
Hamiltonian parameters for STO. 

Our work emphasizes the need for additional exper- 
imental data on the electronic band structure of Per- 
ovskites. Even for STO, by far the most studied d° 
Perovskite, existing experimental data is insufficient to 
uniquely determine the values of band parameters that 
will, for example, control the character of the two- 
dimensional electron systems formed by ^-doping. 

In the past few years much effort has been devoted to 
fabricating oxide-based hetro-structures. Our model for 
the electronic structure of the bulk is a first step towards 
modeling these complex systems. 

This work was supported by Welch Foundation Grant 
F1473. We acknowledge helpful conversations with 
Jim Allen, Young Jun Chang, Harold Hwang, Worawat 
Meevasana, and Susanne Stemmer. 
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FIG. 8: (Color online) Dependence of magnetic oscillation 
frequency and cyclotron mass on A S o- Here L/M = 1/8, 
N = and A T = 0. (a,b): Scaled SdH frequency for magnetic 
fields along [001], [110], and [100] as a function of A S o- (c,d): 
Cyclotron mass as a function of A S o- 



electronic density. Similar graphs are given in Figj8]for 
a scenario in which A T <C A so . The different trends of 
F and m* as a function of density clearly distinguishes 
between the A T ^> A so scenario and its opposite coun- 
terpart. 



V. SUMMARY 

d° Perovskites have played a centeral role in various 
areas of solid state physics and are now emerging as im- 
portant building blocks for oxide-based hetro-structures. 



Appendix A: k • p Hamiltonian for the tetragonal 
phase 



The momentum dependent part of the effective hamil- 
tonian h for a d° perovskite is given by Eq.(§. In 
this appendix we use group theory methods to ex- 
press h in terms of a small number of phenomenological 
parameters^. 

The calculation of h involves the evaluation of matrix 
elements of the form (0|k • p^-). Here ipj is a basis func- 
tion of the tig manifold, p is the momentum operator, 
and is a state outside of the ti g manifold. At the cubic 
to tetragonal phase transition the symmetry at the zone 
center reduces from Oh to D^. Correspondingly, at the 
phase transition the three ^-'s change their transforma- 
tion properties — > T± + . The three components 
of the momentum operator p, that transform as a single 
irrep (IY 5 ) in the cubic phase split: p x ->Py £ whereas 

p z e 17. 

The values of the matrix elements vary smoothly across 
the structural transition. To emphasize the relation be- 
tween the two symmetries we label the matrix elements 
in the tetragonal phase with a subscript that corresponds 
to the irrep of <j) in the cubic phase and a superscript 
that corresponds to its irrep in the tetragonal phase . 
For example, is associated with a basis function that 
evolved from in the cubic phase to T 5 in the tetrag- 
onal one. 



We first consider the TJ band. The intermediate states 



are 



r 3 "- 



(Al) 



9 



Denoting 



(^5x\(Px^Py)\^t ) — 1^(0^1) 

<r 5 - \( Px , Py )\r+) = | B ^ 
(r 3 -|p z |r 4 + ) = { 



o) 



Bf 2 



we find that 



h zz =M B (kl + k 2 y ) + C B k 2 z 



(A2) 



(A3) 



where the two real phenomenological parameters are 
given by 



Mi 



1 

2m ' m 



n 2 ^ 



\Bl, 



5 |2 



and 



1 

2m ' rn 



\Bt, 



3 |2 



, er -E^t(0)-E r J(0) 



(A4) 



(A5) 



We now turn to the band. The intermediate states 



are 



>r 5 + 



(A6) 



Following similar steps to those taken above we obtain ex- 
pression (10). The k-dependent Hamiltonian h depends 



on six real parameters and a single complex parameter 
Mbc- The k • p Hamiltonian in the cubic phase can easily 
be obtained from its tetragonal counterpart by disregard- 
ing the subscripts of the phenomenological parameters; 
for example by associating with Cb and Cc a single pa- 
rameter C. 
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